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ABSTRACT 

Extreme Ultraviolet Explorer Deep Survey observations of cool stars (spectral type F to M) have been used 
to investigate the distribution of coronal flare rates in energy and its relation to activity indicators and rotation 
parameters. Cumulative and differential flare rate distributions were constructed and fitted with different methods. 
Power laws are found to approximately describe the distributions. A trend toward flatter distributions for later- 
type stars is suggested in our sample. Assuming that the power laws continue below the detection limit, we 
have estimated that the superposition of flares with radiated energies of about 10 29 - 10 31 ergs could explain the 
observed radiative power loss of these coronae, while the detected flares are contributing only w 10 %. While the 
power-law index is not correlated with rotation parameters (rotation period, projected rotational velocity, Rossby 
number) and only marginally with the X-ray luminosity, the flare occurrence rate is correlated with all of them. 
The occurrence rate of flares with energies larger than 10 32 ergs is found to be proportional to the average total 
stellar X-ray luminosity. Thus, energetic flares occur more often in X-ray bright stars than in X-ray faint stars. 
The normalized occurrence rate of flares with energies larger than 10 32 ergs increases with increasing Lx/iboi 
and stays constant for saturated stars. A similar saturation is found below a critical Rossby number. The findings 
are discussed in terms of simple statistical flare models in an attempt to explain the previously observed trend for 
higher average coronal temperatures in more active stars. It is concluded that flares can contribute a significant 
amount of energy to coronal heating in active stars. 

Subject headings: stars: activity — stars: coronae — stars: flare — stars: late-type — stars: rotation — X-rays: 
stars 



1. INTRODUCTION 

Stellar activity of "normal" stars has been explored exten- 
sively in the X-ray regime for more than two decades (e.g., Va- 
iana et al. 1981; Linsky 1985). The chromospheres and coro- 
nae of some late-F to M main-sequence stars have been found 
to show enhanced magnetic activity. The latter is underlined by 
enhanced activity indicators such as the X-ray luminosity Lx, 
■ its ratio to the bolometric luminosity Lx/Lboi, the presence of 
flares in the optical U band (and in other wavelength regions), 
flux variations in chromospheric lines, spots on the stellar sur- 
face, etc. A dynamo mechanism is thought to be the primary 
cause for stellar and solar activity as suggested by the empirical 
relation between rotation and activity (for a recent review, see 
Simon 2000). Rotation parameters (such as the rotation period, 
the Rossby number, or the angular velocity) are therefore prime 
parameters that determine the stellar activity level. 

Stellar rotation was proposed to determine the level of activ- 
ity of solar-type stars by Kraft (1967). Quantitative relation- 
ships between activity indicators and rotation parameters (e.g., 
Skumanich 1972; Pallavicini et al. 1981; Walter 1982; Noyes 
et al. 1984; Randich et al. 1996) provide information on the 
physical origin of stellar activity. In a study of X-ray emission 
from stars, Pallavicini et al. (1981) found that X-ray luminosi- 
ties of late-type stars are dependent on the projected rotational 
velocity but are independent of bolometric luminosity. Walter 
& Bowyer (1981) and Walter (1981) presented observational 
evidence that RS CVn systems and G-type stars show a quies- 
cent Lx/^boi ratio proportional to their angular velocity. How- 
ever, in the same series of papers, Walter (1982) proposed that 
the simplistic view of a power-law dependence should be re- 



placed by either a broken power law or by an exponential de- 
pendence of Lx/L\, \ on the angular velocity. This led to the 
concept of saturation of stellar activity (Vilhu 1984; Vilhu & 
Walter 1987) at high rotation rate. Noyes et al. (1984) sug- 
gested that the Rossby number Rq (defined as the ratio P/t c 
of the rotation period P and the convective turnover time r c ) 
mainly determines the surface magnetic activity in lower main- 
sequence stars. Stepieh (1993) showed that, for main-sequence 
late-type stars, Rq correlates with activity indicators better than 
P does. 

Flares are direct evidence of magnetic activity in stellar at- 
mospheres. They are also in the center of the debate on the 
origin of coronal heating. Although several possible heating 
mechanisms have been identified (e.g., Ionson 1985; Narain 
& Ulmschneider 1990; Zirker 1993; Haisch & Schmitt 1996), 
there is increasing evidence that flares act as heating agents of 
the outer atmospheric layers of stars. A correlation between 
the apparently non-flaring ("quiescent") coronal X-ray lumi- 
nosity Lx and the stellar time-averaged f/-band flare luminosity 
(Doyle & Butler 1985; Skumanich 1985) suggests that flares 
can release a sufficient amount of energy to produce the sub- 
sequently observed quiescent coronal emission. Robinson et 
al. (1995, 1999) found evidence for numerous transition region 
(TR) flares in CN Leo and YZ CMi. Further evidence for dy- 
namic heating has been found in broadened TR emission line 
profiles (Linsky & Wood 1994; Wood et al. 1996) that were 
interpreted in terms of a large number of explosive events. 

Stellar X-ray and EUV flares have been found to be dis- 
tributed in energy according to a power law (see Collura, 
Pasquini, & Schmitt 1988; Audard, Gudel, & Guinan 1999; 
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Osten & Brown 1999), similar to solar flares (e.g, Crosby, As- 
chwanden, & Dennis 1993). One finds, for the flare rate dN 
within the energy interval [E,E + dE], 

dN = kl E-«. 



dE 

The cumulative distribution for a > 1 is defined by 
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where the normalization factors k\ and k 2 = k\j (a— 1) are con- 
stants. For a > 2, an extrapolation to flare energies below the 
instrumental detection limit could be sufficient to produce a ra- 
diated power equivalent to the X-ray luminosity Lx of the qui- 
escent corona, 
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where £ m „ is the energy of the most energetic relevant flare. 
Small values of the minimum flare energy E m [ n then can lead 
to an arbitrarily large radiated power. Thus, as pointed out by 
Hudson (1991), it is crucial to investigate whether the solar (and 
stellar) flare rate distributions in energy steepen at lower ener- 
gies. Parker (1988) suggested that the heating of the quiescent 
solar corona could be explained by "microflares". In the solar 
context, several estimates of the power-law index have recently 
been given. For "normal" flares, a ~ 1.5—1.8 (Crosby et al. 
1993); it takes different values for smaller flare energies, from 
a ;=y 1 .6 for small active-region transient brightenings (Shimizu 
1995) to a = 2.3-2.6 for small events in the quiet solar corona 
(Krucker & Benz 1998), and other "intermediate" values (e.g., 
Porter, Fontenla, & Simnett 1995; Aschwanden et al. 2000; Par- 
nell & Jupp 2000). Stellar studies on flare occurrence distribu- 
tions in X-rays are rare, probably due to the paucity of stellar 
flare statistics. Using EXOSAT data, Collura et al. (1988) found 
a power-law index a of 1.52 for soft X-ray flares on M dwarfs. 
Osten & Brown (1999) reported a = 1.6 for flares in RS CVn 
systems observed with EUVE. On the other hand, Audard et al. 
(1999) found a power-law index a w 2.2 ± 0.2 for two young 
active solar analogs. 

In this paper, we present a follow-up study of Audard et 
al. (1999) for EUVE observations of active late-type main- 
sequence stars. Together with an investigation of coronal heat- 
ing by flares, a more general picture of the relation between 
flares and stellar activity in general is developed. Section 2 
presents the method for the data selection and reduction, sec- 
tion 3 explains the construction of cumulative and differen- 
tial flare rate distributions in energy. Section 4 provides the 
methodology for fitting the distributions, while section 5 ex- 
plores quantitatively the correlations between various physical 
parameters. Finally, section 6 gives a discussion of the results, 
together with conclusions. 

2. DATA SELECTION AND REDUCTION 

We use data from the Extreme Ultraviolet Explorer {EUVE, 
e.g. Malina & Bowyer 1991) to study the contribution of flares 

4 For 47 Cas, the X-ray emitter is the probable, optically hidden G0-5 V companion (Gudel et al. 2000), therefore we have set B-V = 0.62. We have estimated the 
equatorial velocity from the rotation period of the X-ray bright source, from the bolometric luminosity and an effective temperature of r eff 5900 K. 

5 IRAF is distributed by the National Optical Astronomy Observatories (NOAO). STScI and NOAO are operated by the Association of Universities for Research 
in Astronomy, Inc. 



to the observable EUV and X-ray emission from stellar coro- 
nae. In order to identify a sufficient number of flares in the 
EUVE Deep Survey (DS) light curves, data sets with more than 
5 days of monitoring, or with a significant number of flares 
(more than ten flares identified by eye) were selected. Active 
coronal sources were our prime choice, as these stars often 
show several distinct stochastic events. We have focused our 
analysis on young, active stars that do not display rotationally 
modulated light curves and that can be considered as single X- 
ray sources. Some stars in the sample are detected or known 
binary systems, in which only one of the components is be- 
lieved to contribute significantly to the EUV and X-ray emitted 
radiation. If a star was observed several times (more than a 
few days apart), we considered the different data sets as orig- 
inating from different coronal sources, as the activity level of 
these stars is usually not identical at two different epochs. We 
carefully checked the DS data and rejected data that presented 
evident problems, such as "ghost" images in the DS remapped 
event files, or incursions into the DS "dead spot". The final list 
contains 12 stellar sources (1 F-type, 4 G-type, 2 K-type, and 
5 M-type coronal sources). We do not claim our sample to be 
complete in any sense. However, this sample is representative 
of the content of the magnetically active cool main-sequence 
stellar population in the EUVE archive. Table 1 gives the name 
of the stellar source (Col. 1), its spectral type (Col. 2), its dis- 
tance d in parsecs from Hipparcos (Perryman et al. 1997; ex- 
cept for AD Leo and CN Leo which are from Gliese & Jahreiss 
1991; Col. 3), the rotation period P in days together with its 
reference (Cols. 4 & 5), the projected rotational velocity 4 and 
its reference (Cols. 6 & 7), the color index B-V 4 and the visual 
magnitude V from Hipparcos (Perryman et al. 1997; except for 
AD Leo and CN Leo for which data were retrieved from Sim- 
bad; Cols. 8 & 9), the mean DS count rate (Col. 10), the derived 
(see below) EUV+X-ray (hereafter "coronal") luminosity in the 
0.01-10 keV energy range (Col. 1 1), and the EUVE observing 
window (Col. 12). 

We have made extensive use of the data from the EUVE 
archive located at the Multimission Archive at the Space Tele- 
scope Science Institute (STScI) 5 . DS Remapped Archive QPOE 
files were rebuilt using the euvl . 8 package within IRAF 5 . 
Light curves (Fig. 1) were created using a DS background 
region ten times larger than the source region area. Event 
lists were extracted from the source region for further analy- 
sis. Thanks to the sufficiently large count rates of our sources, 
the contribution of the DS background was very low, and could 
be neglected (after a check for its constancy). Also, with our 
analysis method, flare-only count rates (count rates above the 
"quiescent" level) were required; therefore the small contribu- 
tion of the background was eliminated in any case. Event and 
Good Time Interval (GTI) files were then read and processed 
with a flare identification code. Audard et al. (1999) explain in 
detail the procedure applied to identify flares in the DS event 
files. In brief, the method, adapted from Robinson et al. (1995), 
performs a statistical identification of flares. It assigns occur- 
rence probabilities to light-curve bins. Note that several time 
bin lengths and time origins for the binning are used so that the 
identification of flares is not dependent on the choice of these 
parameters. Note also that, due to gaps between GTIs, the ef- 
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fective exposure of a bin had to be taken into account. We refer 
to Robinson et al. (1995) and Audard et al. (1999) for more 
details. 

Figure 1 (upper panels) shows background-subtracted EUVE 
DS light curves 6 for some data sets. Only for plotting purposes, 
the data have been binned to one bin per orbit (P or b = 96 min). 
Note that for our data analysis, we have not restricted ourselves 
to the above bin size: bin durations from 1/5 to twice the orbital 
period have been used (see Fig. 1). The lower panels show the 
corresponding "significance plots", which give the probability 
for the presence of quiescent bins as a function of time (x-axis) 
and bin size (y-axis). The flare significance increases from light 
gray to black. Physical parameters (start time, end time and to- 
tal duration) were determined from smoothed high-resolution 
light curves, using a Gaussian fit to the flares above a smooth 
lower envelope characterizing the quiescent contribution. Thus 
the start and end times of a flare were defined as the times sep- 
arated from the maximum by 2 a, where a is the standard de- 
viation of the Gaussian function. Within this interval, we cal- 
culated a mean flare count rate above the quiescent emission by 
subtracting the mean background levels just before and after the 
flare, and multiplied it by the flare total duration to derive the 
total number of flare counts C. We then used a constant count- 
to-energy conversion factor (/ = 1.06 x 10 27 ergs ct" 1 pc~ 2 ) to- 
gether with the source distance d to derive the total energy E 
radiated in the EUV and X-rays, 



E = Cxfx (4ird 2 



(4) 



We derived the conversion factor / from mean DS count 
rates of archival EUVE cool-star data sets and published X- 
ray luminosities (Pallavicini et al. 1988; van den Oord, Mewe, 
& Brinkman 1988; Pallavicini, Tagliaferri, & Stella 1990; 
Dempsey et al. 1993a,b; Schmitt, Fleming, & Giampapa 1995; 
Monsignori Fossi et al. 1996; Dempsey et al. 1997; Tagliaferri 
et al. 1997; Audard et al. 1999; Hunsch et al. 1999; Sciortino et 
al. 1999). We corrected the published X-ray luminosities to the 
new, Hipparcos-derived distances (Perryman et al. 1997). Then, 
using a typical model for young active stars (2-temperature col- 
lisional ionization equilibrium MEKAL model [T\ = 0.6 keV, 
T2 = 2.0 keV] with the iron abundance Fe = 0.3 times the so- 
lar photospheric value), we estimated factors to apply to the 
published luminosities in order to convert them to 0.01-10 keV 
luminosities. We finally defined the conversion factor / as 
the mean ratio between the observed fluxes Lx/(4-ird 2 ) and the 
mean DS count rates /j,. Note that for our targets, the corona 
radiates mostly in the X-ray band (w 0.1-5 keV) rather than in 
the EUV band. 

3. FLARE OCCURRENCE RATE DISTRIBUTIONS 

Cumulative flare rate distributions in energy were con- 
structed for each source (Fig. 2). Similarly to Audard et al. 
(1999), we applied a correction to the effective rate of identi- 
fied flares. For each cumulative distribution, the flare rate at 
the energy of the second-largest flare was corrected by a factor 
Aotai /(Aotai -Aargest), where D tota i is the total observing time 
span, and Aargest is the total duration of the largest flare. Anal- 
ogously, the correction for the third-largest flare took into ac- 
count the total durations of the largest and second-largest flares, 
and so on. This correction was necessary since usually more 
than 50 % of the DS light curves were occupied by identified 
flares, i.e. it is common for flares to overlap in time. 

6 With task qpbin of euvl . 8, the last bin of a light-curve plot is generally 



For each cumulative distribution, we have two series of pa- 
rameters, namely the flare energy (£,) and the flare rate at this 
energy (si), with the indices running from (largest-energy 
flare, so = 1 /Aotai by definition) to M (lowest-energy flare). The 
cumulative occurrence rates N(> Ef), i.e., the rate of flares per 
day with energies exceeding E it were defined as 

N(> Eq) = -J— = s , (5a) 



Aotal 
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We then constructed differential distributions. For the energy 
interval [E i} E i+ i], we defined differential flare occurrence rates 
(«() as 



Ei-E, 



:0,...,M-1. 



(6) 



We calculated uncertainties for the flare occurrence rates per 
unit energy; we assumed that each number of flares (s- = Sj x 
Aotai) with energy Ej has an uncertainty As'j estimated from 
a Poisson distribution. This approximation accounts for the 
larger relative uncertainty (As'j / Sj) of the number of detected 
flares at higher flare energies than at lower flare energies. How- 
ever, instead of setting As'j = we have used the approxima- 
tions proposed by Gehrels (1986), who showed that for small k 
following a Poisson distribution, the 1 a upper error bar can be 

1 /2 

approximated with 1 + (k+ 3 /4) , while the the 1 a lower er- 
ror bar is still approximated by the usual definition y/k. There- 
fore, we have defined As'j as the geometrical mean of the upper 
and lower error bars: 

1/2 



a,; = ^x(i +v /^ 



(7) 

It follows that each differential occurrence rate n, has an uncer- 
tainty An, equal to 

A^/Aotal 



An: = ■ 



Ei-E, 



(8) 



4. FITS TO THE DISTRIBUTIONS 

4. 1 . Cumulative Distributions 

The power-law fitting procedure to the cumulative flare oc- 
currence rate distributions is adapted from Crawford, Jauncey, 
& Murdoch (1970). This method is based on a maximum- 
likelihood (ML) derivation of the best-fit power-law index a. 
The best-fit normalization factors ki (see eq. [2b]) were then 
computed; using equation (3b), the minimum flare energy E m { n 
required for the power law to explain the mean observed radia- 
tive energy loss (Lx) was calculated, assuming that the cumula- 
tive flare occurrence rate distribution in radiated energy follow 
the same power law below the flare energy detection limit: 

l/(2-a) 



h*( 2ll \+e- 
k 2 \a-l ' max 



(9) 



where the coronal luminosity Lx was estimated from \i, the 
mean DS count rate, 

L X = fXXfX (4Trd 2 ). (10) 

Columns 2 and 3 of Table 2 give the power-law indices a of the 
cumulative distributions, and the corresponding minimum flare 
energies E m [ n . The best-fit power-law indices a derived from 
simple linear fits (x 2 method) in the logA^(> E)-logE plane 
have been added for comparison in Column 4. 

omitted (D. J. Christian 1999, priv. comm.) 
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We find a possible trend for decreasing power-law indices 
with increasing color indices; sources of spectral type F or G 
tend to show indices above the critical value of 2, although 
a < 2 is acceptable within the 1 a confidence range (except for 
the F star HD 2726). On the other hand, K and M stars show 
various power-law indices, with a > 2 being usually marginally 
acceptable, although some individual sources show best-fit val- 
ues above 2. The minimum flare energies E m { n can be asso- 
ciated with relatively small stellar flare energies. In the so- 
lar context, however, they correspond to medium-to-large flares 
(E w 10 29 -10 31 ergs). 

4.2. Differential Distributions 

Our cumulative distributions do not account for uncertainties 
in the flare occurrence rates. Differential distributions allow us 
to avoid this effect, and they include uncertainties in a natural 
way (see section 3). Therefore, we propose to use this different 
approach in order to compare the results. The differential distri- 
butions were transformed into FITS files and were read into the 
XSPEC 10.00 software (Arnaud 1996). Due to the characteris- 
tics of XSPEC, energy bins were created in which each bin SE, 
is defined as the interval between two consecutive flare ener- 
gies (namely 6E t = [£,+!, £,], = 0, . . . ,M— 1). Finally, a power- 
law fit (implemented in the software, using the y 2 minimization 
method) was performed for each distribution. To estimate the 
uncertainties derived for the index a, confidence ranges for a 
single parameter were calculated by varying the index a and 
fitting the distribution until the deviation of \ 2 from its best-fit 
value reached dx 2 = 1 00. The power-law indices and their cor- 
responding confidence ranges can be found in Col. 5 (Table 2). 
Note that the small "signal-to-noise" of the distributions did not 
allow us to better determine the confidence ranges. Relative 
uncertainties of the n, values were usually larger than 50 %, 
reaching about 150 % at most. We also note that the values 
of a derived from fits to the cumulative distributions are simi- 
lar to those derived from fits to differential distributions; con- 
fidence ranges for the second method are, however, larger and 
originate from the inclusion of uncertainties in the flare rates, 
together with the small number of detected flares. A strong 
support for this statement comes from the confidence range de- 
rived for k Cet 1994. Only four energy bins were used, leading 
to large confidence ranges and an unconstrained upper limit for 
a. 

4.2.1. Combined Data Sets 

A possible dependence of a with the stellar spectral type has 
been mentioned above. To test this trend further and also to 
obtain tighter results for our power-law fits, we have performed 
simultaneous fits to the differential distributions within XSPEC. 
In brief, all data sets belonging to a spectral type (we combined 
the only F star with the G-type sources) were fitted simultane- 
ously with power laws of identical index a and one normal- 
ization factor for each data set. Note that, with this procedure, 
the confidence range for a is better determined than in the case 
of individual distribution fits. The implicit hypothesis of this 
procedure assumes that the coronae of stellar sources within a 
given spectral class behave similarly, without any influence by 
age, rotation period, projected stellar velocity, etc. Columns 6 
& 7 of Table 2 show the result of the simultaneous fits, together 
with 68.3 % confidence ranges for a single parameter. Again, 

7 logL bo i = logLboi.Q +O.4(M bo i.0 -M bo i), where M bo i = V -(5\ogd-5) + BC. 
Schmidt-Kaler(1982). 



we find a trend for lower indices at later spectral types, although 
the significance is marginal at best. 

5. CORRELATIONS WITH PHYSICAL PARAMETERS 

We have explored correlations of the best-fit power-law in- 
dices a and occurrence rates of flares showing energies larger 
than a typical energy observed in our data (L c = 10 32 ergs) with 
rotation parameters (rotation period P, projected rotational ve- 
locity vsin;', Rossby number Rq) and activity indicators (coro- 
nal luminosity L x , and its ratio L x /Lboi to the bolometric lu- 
minosity). Two nonparametric rank-correlation tests (Spear- 
man's rs and Kendall's r; Press et al. 1992) were used. These 
robust tests allow us to calculate correlation coefficients and to 
obtain a two-sided (correlation or anticorrelation) significance 
for the absence of correlation. Thus, a low correlation coeffi- 
cient (rs or t) can be associated with a high probability that the 
sample is not correlated. Note that Kendall's r is more non- 
parametric than Spearman's rs because it uses only the relative 
ordering of ranks instead of the numerical difference between 
ranks (Press et al. 1992). No uncertainty was included in the 
tests. In Table 3, we give the rank-correlation coefficients to- 
gether with their two-sided significances in parentheses. We 
have defined two data groups for the tests. The first group (here- 
after DG1) corresponds to power-law indices or flare rates de- 
rived from ML fits to the cumulative distributions, while data 
group 2 (hereafter DG2) corresponds to indices or flare rates 
derived from fits to the differential distributions. 

5.1. Correlations of the Power-Law Index a 

5.1.1. Coronal Luminosity Lx 

Coronal luminosities (see eq. [10]) were used to test their 
correlation with a (Fig. 3). For each data group, the signifi- 
cance levels (Table 3) are usually smaller than 5 %, the highest 
level reaching 10 %. Therefore, placing a limit of 5 % to the 
two-sided significance level, the correlation between a and Lx 
is marginally significant at best. Such a correlation can be ex- 
plained as follows. For saturated stars (L x ~ 10~ 3 Lb o i), the X- 
ray luminosity should decline for stars with spectral type from 
F to M (hence increasing color index B—V) because of the de- 
crease of their bolometric luminosity. In section 4.2.1, we have 
found a suggestion that the power-law index a is weakly corre- 
lated with the stellar spectral type. Therefore, a weak correla- 
tion of Lx with a can be expected. Note that due to the scatter in 
the Lx/Lboi ratio in our sample, this can lead to a scatter in the 
dependence of L x on the spectral type, hence in the dependence 
of L x on the index a. 

5.1.2. Ratio Lx/L^oi 

We have tested the correlation between L x /Lb i and the 
power-law index. The bolometric luminosities were calculated 7 
from parameters in Table 1 and corresponding bolometric cor- 
rections. Two-sided significance levels (Table 3) show that the 
correlation between L x /Lt, i and a is not significant. 

5.1.3. Rotation Period P 

For stars with known rotation periods (Table 1), we have 
tested the correlation between the power-law index a and P. 
Note that for n Cet, for which there are two data sets, we have 

Here, Mb o i.0 = 4.64 mag, Lb o i.0 = 3.85 X 10 33 erg s" 1 and BC is taken from 
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used weighted means of the indices a. Our final sample then 
comprised 7 data points. From the two-sided significance lev- 
els, we can state that no correlation between the power-law in- 
dex a and the rotation period P is present in our sample (Ta- 
ble 3). 

5.1.4. Projected Rotational Velocity vsin; 

Projected rotational velocities from Table 1 were used. We 
used the same procedure as above to calculate the weighted 
mean of a for k Cet; for the tests, the upper limits (GJ 411, 
CN Leo) have been omitted. Our sample then contained 8 data 
points for each data group. The two-sided significances for the 
present correlation tests (Table 3) imply the absence of a signif- 
icant correlation between the power-law index a and the pro- 
jected rotational velocity vsin; in our data sample. 

5.1.5. Rossby Number Rq 

We have used the available periods in Table 1 and have calcu- 
lated the convective turnover times t c (2) (the number 2 refers to 
the ratio of the mixing length to the scale height) from the B— V 
color index and equation (4) of Noyes et al. (1984) in order to 
derive Rq. The nonparametric tests again suggest an insignifi- 
cant correlation between the index a and the Rossby number. 

5.2. Correlations of the Flare Occurrence Rate 

5.2.1. Flare Rate vs. L x 

Figure 4 shows the occurrence rate of flares with energies 
larger than 10 32 ergs (E c ) versus the coronal luminosity L X - 
Quiescent X-ray luminosities (corrected to an energy range 
between 0.01 and 10 keV) taken from the literature (Collier 
Cameron et al. 1988; Pallavicini et al. 1990; Hunsch, Schmitt, 
& Voges 1998; Hunsch et al. 1999) were also used for compar- 
ison (crosses), although they were not taken into account in the 
following tests. The evident correlation between N(> E c ) and 
Lx is confirmed by Spearman's r$ and Kendall's r tests. The 
former has rank-correlation coefficients of 0.95 and 0.94, with 
two-sided significances for deviation from zero of 2 x 10~ 6 and 
4 x 10" 6 for DG1 and DG2, respectively. The latter test has 
coefficients of 0.85 and 0.82, while significances are 1 x 10 -4 
and 2 x 10~ 4 . This implies that a correlation between Lx and 
the flare occurrence rate is highly significant for our data sam- 
ple (Figure 4). The linear best-fit in the log-log plane for DG1 
is logM> E c ) = (-26.7 ± 2.9) + (0.95 ± 0. 10) logL x (number of 
flares per day), while the best-fit for DG2 was logM> E c ) = 
(-25.5 ± 2.8) + (0.90 ± 0.10)logL x , hence the relation between 
the flare rate and the luminosity is compatible with propor- 
tionality. Note that, as the correlation between a and Lx was 
marginal at best (section 5. 1 . 1), we can safely state that propor- 
tionality exists between L x and the normalization factor Ic2 (and 
hence ki). 

5.2.2. Normalized Flare Rate vs. L x /L\, i 

The canonical saturation limit for stars with different spec- 
tral types has been found to appear at different v sin ;' (e.g., Cail- 
lault & Helfand 1985; Stauffer et al. 1994; Randich et al. 1996; 
Stauffer et al. 1997), and therefore, based on the relation of 
Pallavicini et al. (1981) between vsin; and L x for unsaturated 
stars, at different X-ray luminosities. Figure 5 shows the coro- 
nal luminosity L x against the bolometric luminosity Lboi- It em- 
phasizes the different loci of our coronal sources with respect to 

8 a = (wi«i +W2Ct2) I (w[ +W2), where w,- = 1/ (<T;,uCi,i), and ct,-, u and try are 



saturation. Note that L x / Lboi ratios range from « 1 0~ 5 to « 1 0" 3 . 
For our correlation tests, we have normalized, for each source, 
the occurrence rates of flares with energies larger than 10 32 ergs 
with the occurrence rate N(> E c ) sat at the saturation turn-on for 
the given spectral type (Lx.sat = 10~ 3 Lb o i) derived from the best- 
fits to N(> E c ) vs. Lx in the previous section. Thus, we are able 
to check whether or not the normalized flare occurrence rate 
stays constant at unity at activity saturation. Figure 6 suggests 
that it does and that the flare rate saturates at activity saturation. 
In the following sections, we will suggest that this effect is not 
biased by the absence of stars "beyond" saturation. Note a few 
discrepant features, such as for 47 Cas (point 2). Its point in 
DG2 is about 1/2 dex higher than in DG1. Since its index is 
larger (a ~ 2.6) for DG2 than for DG1 (a w 2.2), and since its 
minimum observed flare energy is about 10 33 ergs, it follows 
that N(> E c ) is larger for DG2 than for DG1. Similarly, CN 
Leo (points 11 & 12) has two different a indices for the 1994 
and 1995 observations (a ss 2.2 and 1.5). This large discrep- 
ancy (due to the flat low-energy end of the 1995 distributions) 
induces different normalized flare occurrence rates. 

5.2.3. Normalized Flare Rate vs. Normalized vsin; 

As before, it was necessary to normalize the flare occur- 
rence rates. In section 5.2.2, we have shown that, for our 
sample, the normalized flare occurrence rate does not show 
a trend to increase at saturation. However, our sample con- 
tains stars "beyond" saturation, i.e., stars that appear saturated 
(Lx/Lboi ~ 10" 3 ) but that rotate faster than a star at the onset of 
saturation. Therefore, we have normalized v sin ;' with (vsin;') sat , 
where the latter were obtained from the Pallavicini et al. (1981) 
relation (L x = 1.4 x 10 27 [vsin;'] 19 ) at the saturation turn-on for 
each stellar spectral type. Stars at the saturation level thus have 
normalized velocities around 1, while those beyond that level 
have values significantly higher than 1 (Fig. 7). In our sample, 
two stars show high values. The first is the bright Kl dwarf 
AB Dor, and the second is the M6 dwarf CN Leo, although its 
projected rotational velocity is an upper limit. Hence, AB Dor 
is the only star that supports the suggestion of constant normal- 
ized flare rates at vsin; saturation. However, in the following 
section, we will show that the result is also supported by a cor- 
relation with the Rossby number Rq, which is less dependent on 
spectral type. We have performed a fit to the data, assuming a 
saturation function of the type 

N(>E C ) / 1 vsin; \ 

JV(>E c ) sat 6XP ^ C(vsim)saJ' 

with ( as the single fit parameter. We derived £ = 0.56 for 
DG1, and C = 0.33 for DG2. For the fits, we have used the 
estimated equatorial velocity for 47 Cas, the logarithmically av- 
eraged flare rate for n Cet, and we have not taken into account 
the upper limits for CN Leo and GJ 41 1 . 

5.2.4. Normalized Flare Rate vs. Rq 

We have further tested the saturation of the flare rate using 
Rq. Compared to the normalized projected rotational velocity, 
the Rossby number does not contain the uncertainty due to the 
projection angle ;'. Furthermore, it does not need to be normal- 
ized as it already corresponds to a normalized rotation period. 

upper and lower error bars, respectively. For power-law indices of DG1, we have 
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However, the rotation period is not available for each star of 
our sample. Figure 8 (upper and middle panels) suggests a sat- 
uration effect, although again only AB Dor supports it. Also, 
the N(> E c )/N(> £ , c )sat vs. Ro plots (upper and middle panels) 
are very similar to the well-known activity saturation relation 
with the Rossby number (lower panel of Figure 8). Similarly 
to the normalized flare rate saturation, the luminosity satura- 
tion appears at \ogRo < 1.2. Both the normalized flare rate and 
the luminosity decrease above this limit, as previously found 
for L x /L bol (e.g, Randich et al. 1996; Stauffer et al. 1997). 
Two lines were overlaid for log/? () > 1.2. The solid line cor- 
responds to the best fit to the data of Randich et al. (1996, 
logLx/^boi = -4.4- 1.121og/?oX while the dashed line corre- 
sponds to our best-fit solution, logL x /Lboi = (-4.76 ± 0.30)- 
(1.72±0.42)log/?o, with the ratios Lx/^boi of k Cet being log- 
arithmically averaged. Our sample follows approximately the 
relation of Randich et al. (1996). Hence, together with the up- 
per and middle panels, the lower panel of Figure 8 reinforces 
the suggestion that there is a saturation of the flare rate at the 
activity saturation. 

5.3. Flare Power vs. Lx 

In Figure 9, we have plotted the X-radiated power P F from 
the detected flares as a function of the average luminosity Lx- 
There is an obvious correlation (Spearman: rs = 0.93, P[rs = 
0] = 1.2 x 10" 5 ; Kendall: r = 0.82, P[t = 0] = 2.1 x 10" 4 ) be- 
tween the parameters. This again demonstrates the importance 
of the flare contributions to the observed radiation from active 
coronae. Note that, for our sample and for our flare detection 
threshold, about 10 % of the X-ray luminosity originates from 
detected flares. 

6. DISCUSSION AND CONCLUSIONS 

In the present work we have investigated statistical proper- 
ties of EUV flare events on time scales of days and weeks. 
Although the sensitivity of the EUVE DS instrument is quite 
limited, it provides long time series that are sufficient to draw 
rough conclusions on the statistical flare behavior of active stel- 
lar coronae. 

The first aspect of interest is the distribution of flare ener- 
gies. In the solar case, X-ray flares are distributed in energy ac- 
cording to a power law with a power-law index around 2 (e.g., 
Crosby et al. 1993). Many of the active stars studied here are 
quite different from the present-day Sun, with distinct coronal 
behavior. For example, high-energy particles are continuously 
present in their coronae as inferred from their steady gyrosyn- 
chrotron emission (e.g., Linsky & Gary 1983; Giidel 1994); 
quiescent coronal temperatures reach values of 20 MK which 
are typical on the Sun only in rather strong flares. Individual 
flare energies accessible and observed by EUVE in this study 
are, in several cases, not observed on the Sun at all. Since a 
limited amount of magnetic energy is available in the recon- 
nection zones of the coronal magnetic fields, one would even 
expect some upper threshold to observable flare energies from 
solar-like stars (and therefore to the observed flare X-ray radia- 
tion). 

Our investigation on active main-sequence stars has shown 
that, for the energy range observed, the flare occurrence rate 
distributions in energy can be fitted by power laws. We have 
not found significant evidence for broken power laws that indi- 
cate a threshold energy. The largest amount of radiated energy 
was found to be 10 35 ergs in our flare sample, exceeding the 



X-ray output of very large solar flares by two to three orders of 
magnitude. 

Measuring the value of the power-law indices of the distribu- 
tions is pivotal for assessing the role of flares in coronal heating. 
The quite limited statistics make conclusions somewhat tenta- 
tive, although we emphasize the following. The power-law in- 
dices definitely cluster around a value of 2; they may be slightly 
different for different stellar spectral types (Table 2). The stellar 
distributions are thus broadly equivalent to solar distributions, 
which implies that a) the cause of flare initiation in magneti- 
cally very active stars may be similar to the Sun, and b) that 
the trend continues up to energies at least two orders of magni- 
tude higher than observed on the Sun. We are conversely mo- 
tivated to extrapolate our distributions to lower energies given 
the rather large range over which solar flares follow a power 
law. Caution is in order, however, toward low-energy flares for 
which the distributions may steepen (Krucker & Benz 1998). 

Hudson (1991) argued that, for a power-law index above 2, 
an extrapolation to small flare energies could explain the radi- 
ated power of the solar corona. We derived (eq. [9]) the min- 
imum flare energies E m i n required for the power laws to ex- 
plain our stellar X-ray luminosities (Table 2). For stars with 
a > 2 or just barely below 2, minimum flare energies around 
10 29 - 10 31 ergs were obtained. Such energies correspond to 
intermediate solar flares. Explicit measurements of flare ener- 
gies below our detection thresholds will however be required to 
conclusively estimate their contribution to the overall radiation. 

We have found a trend for a flattening of the flare rate dis- 
tributions in energy toward later spectral types. F and G-type 
stars tend toward power-law indices > 2, while K and M dwarfs 
tend toward indices < 2. If supported by further, more sensitive 
surveys, it suggests that flares play a more dominant role in 
the heating of F and G-type coronae, while they cannot provide 
sufficient energy to explain the observed radiation losses in K 
and M dwarfs. On the other hand, part of this trend could be 
due to the bias introduced by the identification method and the 
length of the GTIs. Although not found in our data, later-type 
(K and M-type) stars may show flares that are typically shorter 
than those of G dwarfs, partly because of the smaller distances 
of the stars that give access to less energetic (and therefore, as 
on the Sun, typically shorter) flares. But then the flare dura- 
tion is smaller than the typical GTI gaps (about 3000 s) so that 
flares that occur between the GTIs remain completely unde- 
tected. This effect can considerably flatten the flare rate distri- 
butions. 

Given that our study is restricted to flares with energies typ- 
ically exceeding 10 31 - 10 32 ergs, it is little suprising that the 
observed flare radiation amounts to only a fraction of the to- 
tal EUV and X-ray losses. We infer an (observed) fraction of 
approximately 10 % relative to the average (quiescent) coronal 
luminosity. This lower limit will undoubtedly increase with bet- 
ter instrument sensitivity. Our study therefore clearly indicates 
that flares provide an important and significant contribution to 
the overall heating of active stellar coronae. 

We have further explored whether statistical flare properties 
are correlated with some physical properties of the stars, such as 
activity indicators and rotation parameters. The power-law in- 
dex a does not correlate with any of the rotation parameters (P, 
vsin/, Rq) nor with the ratio Lx/Lbol- The absence of clear cor- 
relations suggest that the activity phenomena related to flares 
are similar on stars of all activity levels. A marginally signif- 
icant correlation with the coronal luminosity was found. This 
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result is probably related to the trend found for the dependence 
of a on the stellar spectral type. 

On the other hand, the flare occurrence rate above a given 
lower energy threshold is correlated with each of the activ- 
ity indicators and rotation parameters. A single power law 
(N[> E c ] oc L x ) fits the correlation between flare rate and the 
coronal luminosity quite well, indicating that energetic flares 
occur more frequently in X-ray luminous stars than in X-ray 
weak stars. In order to compare the activity levels between 
stars of our sample, we have normalized the flare rate to its 
value that a star adopts at its saturation level. The normalized 
flare occurrence rate increases with increasing activity but stays 
constant for saturated stars. Flare rate saturation underlines the 
close relation between flares and the overall "quiescent" coro- 
nal emission. We now ask more specifically whether the appar- 
ently quiescent X-ray radiation could be related to the derived 
flare distributions. 

In simple terms, we expect a larger magnetic filling factor on 
magnetically more active stars, or more numerous (or larger) 
active regions than on low-activity stars. A higher filling factor 
naturally implies a proportionally higher quiescent X-ray lu- 
minosity and proportionally more numerous flares, so that we 
expect a linear correlation between Lx and the flare rate. But 
what is the nature of the quiescent emission? 

It is known from X-ray observations that the average qui- 
escent coronal temperatures characteristically increase with in- 
creasing activity. The concept of an average coronal temper- 
ature is, however, somewhat problematic. Schrijver, Mewe, 
& Walter (1984) used single-F fits to Einstein/lPC data that 
roughly imply that the total stellar volume emission measure 
in X-rays EM* is proportional to T 3 (as discussed in Jordan & 
Montesinos 1991). Since the radiative cooling function A(T) 
in the range of interest (3 - 30 MK) scales approximately like 
T'^ with <j) « 0.3 (Kaastra, Mewe, & Nieuwenhuijzen 1996), 
we have for L x 

L x ~ EM* A(T) oc EM* T~^ (12) 

and hence Lx oc T ln . Considering that the typical coronal emis- 
sion measure is distributed in temperature, multi-r fits appear 
to better represent average coronal temperatures. Also, to dis- 
entangle functional dependencies from other stellar parameters 
(e.g., the stellar radius), a uniform sample of stars should be 
used. Giidel et al. (1997) derived two-temperature models from 
ROSAT data for a sample of stars that differ only in their activ- 
ity levels but are otherwise analogs to the Sun. For the hotter 
component, they found Lx oc Tj^ 1 6-5 09 , the range of the expo- 
nent illustrating two different spectral models applied. If, how- 
ever, we compute the mean of the two temperatures weighted 
with the corresponding emission measures from their Table 3, 
we find L x oc T 4 8 independent of the spectral model. The two 
more active stars for which ASCA 3-T fits were available in 
Giidel et al. (1997) agree well with this trend (for cooler coro- 
nae, ASCA is not sufficiently sensitive to derive a coronal emis- 
sion measure distribution). Note that the 3-T fits in turn rep- 
resent the derived emission measure distributions in Giidel et 
al. (1997) quite well. We therefore conclude that the emission- 
measure weighted average coronal temperature roughly scales 
as L x oc T 4 - 8 , although with a considerable uncertainty in the 
exponent. 

Hearn (1975) and Jordan et al. (1987) discuss "minimum 
flux" coronae that follow a relation EM* oc T 3 g+ where g* is 
the stellar surface gravity. In a similar way, Rosner, Tucker, 
& Vaiana (1978) find scaling laws for closed static magnetic 



loops that relate external heating, loop pressure, loop-top tem- 
perature, and loop length. The two scaling laws combine to 

y21/6 



Eh oc 



L 2 



(13) 



where E n is the heating rate, T is the (dominating) loop-top 
temperature, and L is the loop length. Schrijver et al. (1984), 
based on these scaling laws, conclude that different families of 
loops must be present on active stars. In either case, explaining 
the locus of the measured (T, EMJ on the empirical relation 
requires an explanation for a specific amount of heating energy 
input to the system. Our investigation suggests that flare events 
contribute significantly to the observed overall emission. If the 
concept of a truly quiescent radiation is abandoned completely 
for magnetically active stars, then we may ask whether flares 
themselves can explain the observed relation between EM* and 
T while at the same time accounting for the proportionality be- 
tween flare rate and L X - We briefly discuss two extreme cases: 

i) Filling factor-related activity. We first assume that flares 
are statistically independent and occur at many different, un- 
related flare sites. The quiescent emission in this case is the 
superposition of all flare light curves. The proportionality be- 
tween flare rate and L x is simply due to more numerous flare 
loops on stars with higher magnetic filling factors (and hence 
higher activity), in which the explosive energy releases build 
up the observed loop emission measure. We ask whether the 
statistical flare distribution determines the average coronal tem- 
perature as well. 

From our flare samples, we have found that the flare dura- 
tion does not obviously depend on the flare energy (as typically 
found for impulsive flares on the Sun; Crosby et al. 1993). We 
therefore first assume one fixed time constant for all flares. For 
example, radiatively decaying flares with a similar evolution 
of coronal densities, abundances and temperatures show sim- 
ilar decay time scales. In this case, the single flare peak lu- 
minosity Cx is proportional to the total radiated flare energy 
E. Further, for solar (and stellar) flares, a rough relation be- 
tween peak emission measure EM and peak temperature T has 
been reported by Feldman, Laming, & Doschek (1995). For 
the interesting range between 5-30 MK, the relation can be 
approximated by 

(14) 



EM = aT h [cm -3 ] 



with a rts 10 cm , b w 5 ± 1, and T measured in K. We ap- 
proximate the flare contribution to the radiative losses by the 
values around flare peak. For the flare radiative losses, we have, 
equivalent to equation (12), 

C x ~ EM A(T) oc EM T~^ (15) 
again with <f> w 0.3 over the temperature range of interest. With 
equations (14) and (15), we obtain 

E oc Lx oc T^. (16) 
To derive a characteristic, emission-measure weighted, time- 
averaged mean coronal temperature T, we average over all flare 
temperatures [T s ,To] by using their peak EM and their occur- 
rence rates as weights: 

L r ° T EM(T) dN IdE dEldT dT 

T = ^= . (17) 

EM(T) dN/dE dE /dT dT 

Here, T s represents the typical temperature of the smallest con- 
tributing flares. From equations (14)— (17), we obtain 



z+ 1 



-(T,/T Q f 



z+l 



i-(T s /T y 



(18) 
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where 

z=(2-a)b-(l-a)(f>. (19) 

Without loss of generality, we assume T s = 1 MK. The up- 
per temperature limit 7b corresponds to the largest flare energy 
Eq typically contributing to the apparently quiescent emission. 
Since the flare rate is given by the power law (eq. [1]), the char- 

1 la. 

acteristic value for Eo scales like fc/ , i.e., with equation (16), 

Tq oc k l J (a(b ~^ . The normalization factor k\ is proportional to 
the overall stellar X-ray luminosity Lx for given a. The steep- 
est dependence of T on k\ is obtained in the limit of small a, 
i.e. z > 0, so that a < (2b-<f>)/(b-<j)), implying a < 2.1. Then, 

for 7b ^> T s , we find T cc 7b, i.e., Lx oc with j3 = a{b — (j)). 
For reasonable values of a = 1.5-2.1, we find f3 w 5.5 — 12, i.e., 
a dependence that is at least somewhat steeper than observed. 

We can repeat the derivation under the assumption of some 
dependence between flare energy E and duration D. Typi- 
cally, larger flares last longer. The function D{T) = D(E[T]) 
should then be used as an additional weighting factor in equa- 
tion (17). Assuming that D cc E x l 2 (e.g., in the case of an 
energy-independent "flare curve shape"), we find E oc C x . In 
that case, equation (18) remains valid, while z = (4 - 2a)b - (3 - 

2a)4> and Lx oc T 13 with (3 = 2a(b-<j>) for "small" a, hence a 
steeper dependence on T than for the first case. 

ii) Loop reheating. In the other extreme case, the flares re- 
peatedly reheat the same coronal plasma in certain "active" loop 
systems (e.g., Kopp & Poletto 1993). In that case, for most 
of the time the coronal loops fulfill a quasi-static approxima- 
tion equivalent to the loop scaling law given by Rosner et al. 
(1978) (see Jakimiec et al. 1992). The situation is equivalent 
to steadily heated loops. The first scaling law of Rosner et al. 
(1978), r max oc (/?L) 1/3 , implies EM oc T 4 and hence L x oc T 3J 
for a given loop length with T in the range of interest. For 
loops exceeding the coronal pressure scale height, Schrijver et 
al. (1984) give an expression EM oc T 5 (with other parameters 
fixed), and hence Lx oc T AJ . In this limit, we attribute the higher 
temperatures of more active stars to higher reheating rates in a 



similar number of loops rather than to a larger number of sta- 
tistically independent heating events in more numerous active 
regions or loops. The cause of an enhanced reheating rate on 
more active stars in a similar number of loop systems remains 
to be explained, however. 

Comparing the observed average coronal temperatures with 
the two extreme results, we see that the observed trend lies in 
the middle between the extreme values. Flare heating is thus a 
viable candidate to explain the trends seen in L x and in T, al- 
though it is not conclusive whether flares occur independently 
or whether they act as reheating agents of coronal loops. It 
appears unlikely that either extreme is appropriate. While a 
larger magnetic filling factor on more active stars undoubtedly 
produces more numerous active regions and thus a higher flare 
rate, the higher coronal filling factor is also likely to lead to 
more numerous reheating events. This is compatible with the 
observed trend between L x and T. 

This study has shown that flares can provide a significant 
amount of energy to heat the coronae of active stars. Although 
the definite answer to which mechanism is responsible for coro- 
nal heating is not yet available, our data sample suggests that 
flares are promising contributors. Better sensitivity together 
with uninterrupted observations (such as those provided by the 
new generation of X-ray satellites XMM-Newton and Chandra) 
will be needed to solve part of the mystery. Our method is 
limited to flares explicitly detectable in the light curves. Al- 
ternative methods are available that model light curves based 
on statistical models. These investigations are the subject of a 
forthcoming paper (Drake et al. 2000). 
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support from the Swiss National Science Foundation, grant 21- 
49343.96, from the Swiss Academy of Sciences and from the 
Swiss Commission for Astronomy. 
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FIG. 1. — Examples of light curves and significance plots. 

FIG. 2. — The cumulative flare occurrence rate distributions for all twelve data sets. 
FIG. 3. — Power-law index a vs. coronal Lx. 

FIG. 4. — Occurrence rate of flares with energies larger than 10 32 ergs vs. coronal luminosity Lx- Data group 1 corresponds to the rates derived from the fits to 
cumulative distributions, while data group 2 corresponds to the rates derived from the fits to differential distributions. The \ 2 linear best-fits in the log— log plane 
are shown as straight lines together with the analytical formulation. Crosses represent Lx values for these sources taken from the literature (see text). 

FIG. 5. — Coronal luminosity Lx vs. bolometric luminosity Lboi- The source identification numbers refer to Table 1. The solid line represents the saturation level 
(Lx/L bo i = 10~ 3 ), while the dotted line is for Lx/L bo i = 10 -4 , and the dashed line for Lx/L bo i = 1CT 5 . 

FIG. 6. — Normalized occurrence rate of flares with energies larger than 10 32 ergs vs. ratio Lx/L bo i of the coronal luminosity and the bolometric luminosity. Data 
groups 1 and 2 are identical to Fig. 4. The numbers are as in Fig. 5. The dash-dotted lines are lines with slope 1. The dotted lines represent the saturation level 
(Lx = 10~ 3 Lboi). Note that, for clarity, the N(> £ c ) of points 4 and 9 have been multiplied and divided by a factor of 1.5, respectively. 

FIG. 7. — Normalized occurrence rate of flares with energies larger than 10 32 ergs vs. normalized projected rotational velocity v sin i/(v sin i)sat =x (spectral-type 
dependent). Data groups 1 and 2 are identical to Fig. 4. "Saturation" fits of the form 1 -exp(-x/Q are plotted dot-dashed. Arrows indicate upper limits. Note that, 
for clarity, the N(> £ c ) of points 4 and 5 have been divided and multiplied by a factor of 1.5, respectively. 

FIG. 8. — (Upper and middle panels) Normalized flare rate vs. Rossby number. Data groups and identifications are as in Fig. 5. (Lower panel) Ratio Lx/L bo i vs. 
the Rossby number for our sample. The solid line corresponds to the fit by Randich et al. (1996), while the dashed line refers to our best fit (see text). 

FIG. 9. — X-radiated power Pf from the detected flares vs. coronal luminosity Lx- The solid line represents proportionality (Pf = Lx), while the dashed line is for 
/*; 0.1L X . 
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Table 1 



Target Selection List. 



Source 


Spectral 


d 


P 


Refi 


v sin i 


Ref 2 


B-V 


V 


A* 


logLx 


EUVE 


Name 


Type 


(pc) 


(d) 




(kms^ 1 ) 




(mag) 


(mag) 


(cts" 1 ) 


(erg s~') 


Observing Window 


HD 2726 


F2 V 


45.07 






13.2 


1 


0.367 


5.67 


0.11 


30.47 


1995/08/09 ->■ 1995/08/16 


47 Cas 


GO-5 V 


33.56 


1.0 


1 


62.1" 




0.620" 


a 


0.14 


30.31 


1997/01/23 ->■ 1997/01/29 


EK Dra 


G1.5 V 


33.94 


2.605 


2 


17.3 


2 


0.626 


7.60 


0.08 


30.09 


1995/12/06 1995/12/13 


k Cet 


G5 V 


9.16 


9.4 


3 


3.9 


3 


0.681 


4.84 


0.08 


28.98 


1994/10/13 1994/10/18 




















0.11 


29.10 


1995/10/06 1995/10/13 


AB Dor 


Kl V 


USA 


0.515 


4 


93.0 


4 


0.830 


6.88 


0.28 


30.05 


1994/11/12 -> 1994/11/17 


e Eri 


K2 V 


3.22 


11.3 


5 


2.0 


3 


0.881 


3.72 


0.30 


28.62 


1995/09/05 -<■ 1995/09/13 


GJ411 


M2 V 


2.55 






<2.9 


5 


1.502 


7.49 


0.02 


27.29 


1995/03/22 -> 1995/04/04 


AD Leo 


M3 V 


4.90 


2.7 


6 


6.2 


5 


1.540 


9.43 


0.28 


28.95 


1996/05/03 -<■ 1996/05/06 


EV Lac 


M4.5 V 


5.05 


4.376 


7 


6.9 


5 


1.540 


10.29 


0.16 


28.74 


1993/09/09 ->■ 1993/09/13 


CN Leo 


M6 V 


2.39 






<2.9 


5 


2.000 


13.54 


0.03 


27.38 


1994/12/16 -> 1994/12/19 




















0.02 


27.27 


1995/01/24 1995/01/30 



a vsim: Estimated equatorial velocity; B-V: Value set to account for the X-ray star's spectral type; V: no data for X-ray emitter; (see text) 



References. — Rotation period (Refi). (1) Giidel, Schmitt, & Benz 1995, (2) Strassmeier, Bartus, & Rodono 1997, (3) Noyes et al. 1984, 

(4) Innis et al. 1988, (5) Baliunas et al. 1983, (6) Spiesman & Hawley 1986, (7) Pettersen, Olah, & Sandmann 1992. Projected rotational 
velocity (Ref 2 ). (1) Groot, Piters, & van Paradijs 1996, (2) Strassmeier & Rice 1998, (3) Fekel 1997, (4) Kiirster, Schmitt, & Cutispoto 1994, 

(5) Delfosse et al. 1998 
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Table 2 

Fits To The Flare Rate Distributions in Energy. 



Cumulative Differential Simultaneous 



Name 


a a 


log(£ min ) d [ergs] 


a h 




a c 


Type 




HD 2726 


2.61 ±0.38 


31.7 (29.7,32.3) 


1.89 e 


2.43 


(1.80,3.93) 


F+G 


2.28 (2.03,2.57) 


47 Cas 


2.19±0.34 


29.7 ( ••• ,31.6) 


1.98 


2.62 


(1.72,5.41) 








EK Dra 


2.08 ±0.34 


30.2 ( • • • , 32.0) 


2.27 


1.78 


(1.26, 2.49) 








k Cet 1994 


2.18±0.89 


27.2 ( ••• ,31.0) 


1.90 


2.55 


(0.31, ••• ) 








k Cet 1995 


2.29±0.51 


29.5 ( ••• ,31.1) 


2.21 


2.45 


(1.65,3.80) 








AB Dor 


1.88 ±0.26 


• • • ( • • • , 28.8) 


1.97 


1.76 


(1.24,2.72) 


K 


1.87 (1.50,2.39) 


e Eri 


2.40±0.81 


29.1 ( ••• ,30.7) 


2.50 


2.38 


(1.06, 4.05) 








GJ411 


1.63 ±0.29 


... ( ... , ... ) 


1.96 


1.57 


(1.08,2.22) 


M 


1.84 (1.63,2.06) 


AD Leo 


2.02 ±0.28 


26.2 ( ••• ,29.8) 


1.85 


1.65 


(1.18,2.35) 








EV Lac 


1.76±0.33 


••• (•••,29.1) 


1.90 


1.75 


(0.98, 3.33) 








CN Leo 1994 


2.21 ±0.30 


29.3 (27.0,29.8) 


1.91 


2.24 


(1.78,3.04) 








CN Leo 1995 


1.46±0.39 f 


... ( ... , ... ) 


2.14 


1.59 


(0.84, 2.56) 









a From an adapted version of Crawford et al. (1970) 
h \ 2 linear fit in the log- log plane 

°X 2 fit within XSPEC with 68 % confidence ranges for a single parameter 

d Minimum energy E m [ n required for the power law to explain the total observed radiative energy loss; limits 
are given in parentheses 

influenced by the largest flare energy; a = 2.43 if removed 
influenced by the flat low-energy end of the distribution 
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Table 3 

Correlation tests for the power-law index a. 



a versus: 




Lx/Lbol 


P 


vsin/ 


Ro 


Spearman test: a 












DG1 


0.50 (0.10) 


-0.28 (0.38) 


0.50 (0.25) 


-0.36 (0.38) 


0.54 (0.22) 


DG2 


0.63 (0.03) 


-0.34 (0.28) 


0.07 (0.88) 


0.05 (0.91) 


0.18(0.70) 


Kendall's r test: a 












DG1 


0.42 (0.05) 


-0.24 (0.27) 


0.33 (0.29) 


-0.29 (0.32) 


0.43 (0.17) 


DG2 


0.45 (0.04) 


-0.33 (0.13) 


0.05 (0.88) 


0.07 (0.80) 


0.14(0.65) 



a Spearman's rank-correlation coefficients rs and Kendall's coefficients r. In paren- 
theses, the two-sided significances of their deviation from zero (=no correlation) 
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FIG. 1 . — Examples of light curves and significance plots. 
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FIG. 2. — Cumulative flare occurrence rate constructed for all twelve data sets. 
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Fig. 3.— 



Power-law index a vs. coronal luminosity Lx. 
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FIG. 4. — Occurrence rates of flares with energies larger than 10 ergs vs. coronal luminosity Lx. Data group 1 corresponds to the rates derived from the fits to 
cumulative distributions, while data group 2 corresponds to the rates derived from the fits to differential distributions. The \ 2 linear best-fits in the log-log plane 
are shown as straight lines together with the analytical formulation. Crosses represent Lx values for these sources taken from the literature (see text). 
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FIG. 5. — Coronal luminosity Lx vs. bolometric luminosity L bo i. The source identification numbers refer to Table 1. The solid line represents the saturation level 
(Lx/Z-boi = 10~ 3 ), while the dotted line is for Lx/L bol = KT 4 , and the dashed line for Lx/L hol = 10~ 5 . 
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FIG. 6. — Normalized occurrence rate of flares with energies larger than 10 32 ergs vs. ratio Lx/Lbol of the coronal luminosity and the bolometric luminosity. Data 
groups 1 and 2 are identical to Fig. 4. The numbers are as in Fig. 5. The dash-dotted lines are lines with slope 1. The dotted lines represent the saturation level 
(Lx = lCT 3 L bo i). Note that, for clarity, the N(> E c ) of points 4 and 9 have been multiplied and divided by a factor of 1.5, respectively. 
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FIG. 7. — Normalized occurrence rate of flares with energies larger than 1() 32 ergs vs. normalized projected rotational velocity v sin j/(v sin i)sat = x (spectral-type 
dependent). Data groups 1 and 2 are identical to Fig. 4. "Saturation" fits of the form 1 -exp(-x/Q are plotted dot-dashed. Arrows indicate upper limits. Note that, 
for clarity, the N(> E c ) of points 4 and 5 have been divided and multiplied by a factor of 1.5, respectively. 
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FIG. 8. — (Upper and middle panels) Normalized flare rate vs. Rossby number. Data groups and identifications are as in Fig. 5. (Lower panel) Ratio Lx/L hol vs. 
the Rossby number for our sample. The solid line corresponds to the fit by Randich et al. (1996), while the dashed line refers to our best fit (see text). 
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